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I present quenched, lattice QCD calculations of the hadronic matrix elements relevant for 






-B°(5) and K" — K'^ mixing and briefly review the status of lattice predictions. 



1 Introduction 

Neutral meson mixing is a rich source of information on the Standard Model (SM). For instance, the 
frequencies with which Bd and Bs mesons oscillate into their anti-particles yield constraints on the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix element Vtd which determines the most poorly known 
side of the unitarity triangle. K^ — K^ mixing, on the other hand, through its measured contribution 
to indirect CP violation in if — > vrvr decays, provides a constraint on the triangle's summit. These 
constraints require quantification of the non-perturbative QCD dynamics which modify the simple, 
underlying quark processes. The uncertainties in this quantification must be reduced to allow for as 
stringent a test of the SM as possible with the triangle's angles soon to be measured at the i?-Factories, 
HERA, the Tevatron and the LHC. Lattice QCD provides a first principle tool which can help achieve 
this goal. In what follows, I present results of lattice calculations performed with C.J. David Lin and 
the UKQCD Collaboration as well as a summary of lattice predictions ^. 



2 RO 



Bd,(s) mixing 



B^ and B^ {Q=d, s) are not eigenstates of the weak hamiltonian and can therefore oscillate into 
one another with a frequency given by the mass difference, Am,g, of the eigenstates oi the full SM 
hamiltonian. In the SM, the dominant contribution to this mass difference is given byt 
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[h''{l-l')q][h,{l-l5)q] 
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Figure 1: Constraints on the Wolfenstein parameters p and r) from Arrid, Arris/ Arriii and e^ Q (for illustration only). 

Xq = {rn^/M^f) and ry^, Sq and C^ are short-distance quantities, calculated perturbatively. Thus, 
\Vtq\ can be determined from a measurement of Arrig once the non-perturbative matrix element 
{Bq\0^^^'^{lJL)\Bq) is quantified ||. This is where the lattice enters. 

While the measurement of Atti^^ provides a direct determination of \Vtd\ 
the ratio 



one may also consider 
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This ratio gives another possible constraint on \Vtd\ since with three generations, \Vts\ — \Vcb\- It 
further has the advantage that many common factors and uncertainties in the evaluation of the matrix 
elements cancel. Measuring Arrr^, however, remains an experimental challenge, as the neutral B_ 
mesons oscillate rapidly: Am^ > 12.4 ps~^ at 95% CL versus Am^ = 0.477(17) ps~^ for B^ mesons 
Nevertheless, even this lower bound on /Arris provides a significant constraint on IV^^^I as can be seen 
in Figure |l|. 
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Traditionally, the matrix element {Bq\Oq '^{p)\Bq) is normalized by its vacuum saturation value: 



BbM 
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imMifi 



(3) 



While one can actually determine the matrix element itself on the lattice, i?-parameters are obtained 
from ratios of correlation functions in which many statistical and systematic uncertainties are expected 
to cancel. Furthermore, the matrix element has mass dimension four and therefore suffers very strongly 
from the uncertainty associated with the determination of the lattice cutoff which is of order 10% in 
present day quenched calculations. As we shall also see later, it is advantageous to get the matrix 
element from an independent determination of Bb^ and /b, and the experimental value of Mb^- 



2.1 A parte on decay constants 

Because the leptonic decay constants of B mesons are required, I briefly digress to comment on their 
values. Marn^ lattice groups have calculated these constants over the years. A recent compilation can 
be found inB where the following summary numbers, which include uncertainties due to quenching, 
are given: 



fBa = 175 ± 35 MeV , Jb. = 200 ± 35 MeV and /b./Zb^ = 1.14 ± 0.08 , 



(4) 



'^ Assuming three-generation unitarity and present day constraints on CKM parameters, |Vtb| = 1 to high accuracy. 
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Table 1: Results for B-meson B-parameters obtained with "relativistic" heavy quarks. (5 is the coupling at which the 

calculations were performed. 13 = oo corresponds to results extrapolated to the continuum limit. ^ is the matching 

scale used. The numbers in italics lare derived from the published results. Running is performed at two-loops using the 

procedure of J. Flynn et al. a which assumes mt — 5 GeV. Bs? is the RG-invariant B-parameter at NLO. 
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in a normalization where fj^ = 131 MeV. While the effects of quenching in (/ss//Bjjl-,appear to be 
small in simulations @, Quenched xPT (Q^PT) indicates that they could be significant 0. 



2.2 AB = 2 matrix elements and B -parameters 

In Table ||, I present our results for Bb^ and Bb^ /Bb^ along with a compilation of results obtained by 
other groups who use, as we do, "relativistic" heavy quarks, as opposed to NRQCD or static quarks. 
(Details of our calculations can be found in the Appendix.) 

Quenching errors for these 5-parameters and 5-parameter ratios have been studied with the help 
of QxPTQ and have been found to be small. Combining this information with the results of Table |l|, 
I give the following estimates: 
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In order to use these results to extract iVj^l from a measurement of Am^, we need to combine 



them with a determination of /b^. Using the estimate given in Section 2A_, I quote: 



fBd^B_ 



'^t 



207(42) MeV 



(6) 



This prediction can be compared to the value obtained from an overconstrained, unitarity-triangle 
fit with /B^y^^ left as a fit parameter 111: /Bd^^^ = 223(13) MeV. This fit incorporates lattice 
predictions for r^d and Bx consistent with the ones given below. Agreement is excellent, indicating 
a general consistency of the SM and lattice calculations. The central value and error bars, of course, 
reflect the choices made by the authors for the various inputs they use. 
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2.3 SU{3) breaking in -^^^ 

There are at least two possible ways of obtaining rg^ from the lattice: 



a) taking the product 



,W 



'MbX ffBsV 
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(7) 



with ifBs/fBa) ^^^ [Bbs/Bb,) determined on the lattice and (MbJMb,) measured experimen- 
tally; 



Table 2: Results for r^d as obtained using methods a) and b) with "relativistic" heavy quarks. 
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b) from a direct determination of the ratio 
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Our results for r^d, together with the results of other groups who use "relativistic" heavy quarks are 
summarized in Table §. Comparison of rj^ at our two values of the lattice spacing (/? = 6.2 and 6.0) 



(a) 



Jb) 



suggests that discretization errors are small. Furthermore, we find that r),J and r),J are compatible, 
though the latter is less accurate and less reliable: its heavy-quark and light-quark-mass dependences 
are stronger and the corresponding extrapolations are less well controlled. 



On the basis of these results and the comments on quenching in Sections 2.1 and 2.2, I quote as 
summary values: 

'Mb, 
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(9) 



3 i^° - K° mixing 

K^ — K^ mixing induces indirect CP violation vn K ^ vrvr decays, quantified by the parameter ex °'- 



eK e-*4 ~ C, CK{lj)BK{lj)A'\'''fi (1 - p)A^r?25o(xi) + Po{xux,, . . .) = (2.280 ± 0.013) x 10 



1-3 



(10) 



with 



(i?°|OA5=2(^)|i^°) = iMlfl X i?x(M) 



and 



Oa5= 



[s-7;.(l-7')rf][s7^1-7')d]. (11) 



This in turn leads to a hyperbolic constraint on the summit (p, ff) of the unitarity triangle, once the 
-B-parameter Bk is determined (see Figure ||). Here, C^ is obtained from well mesured quantities, 
A and A are Wolfenstein parameters and r]2, Ck, Sq and Pq incorporate perturbative, short-distance 
physics {Pq also contains CKM factors). We calculate Bk on the same lattices as the AB = 2 matrix 
elements. 



3.1 Chiral subtractions 

Even though the basic ingredients, such as the operator mixing alluded to in the Appendix, are very 
similar to those used to calculate the AS = 2 matrix elements, the physics here is very different, 
as it is governed by chiral symmetry. In the continuum, Oa5=2 is in the (27, 1) representation of 
SU{3)l X SU{3)r. On the lattice, however, the explicit breaking of chiral symmetry implies the 
following chiral expansion: 



{K^{q)\OAS=2\K'>{p))i,t = aK + pKMl + ^K{p-q) + --- 



(12) 



Table 3: Results for 3}^ (2 GeV). The results in the second half of the table were obtained with discretizations of 
the quark action which maintain a partial or full chiral symmetry, obviating the need for chiral subtractions. 
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where ax and Pk are pm'e lattice artefacts, while ^k{P'Q) and higher-order terms contain the physical 
contributions. In our calculation, where we match onto the continuum at one loop, the artefacts ax 
and fix are proportional to a1 and aa^. The problem is that even though these factors are small, the 
physical contributions are chirally suppressed compared to ax- 

To quantify and subtract the unphysical contributions, we study the chiral behavior of the AS" = 2 
matrix element as a function of M"^ andp-q. At /? = 6.2 we find that artefacts such as ax and j3x are 
small and consistent with zero for all matching scales in the range 1/a — > vr/a. We have checked that 
our results are robust to procedure by normalizing the AS* = 2 matrix element in a variety of ways 
and using different mass and recoil variables for the chiral expansion. The determination of Bx from 
the corresponding physical expansion terms should thus be reliable. At /? = 6.0, the lattice artefacts 
are around 2 standard deviations away from zero and the results are less robust to procedure. Our 
findings, together with results obtained with less improved actions, suggest that discretization errors 
represent an important part of the traditionally observed residual chiral violations. 

3.2 Results for Bx 

We take our [3 = 6.2 result as our best estimate. We run it to 2 GeV at two-loops with n/ = in the 
MS-NDR scheme (small running). Our results, together with those jiif other groups are summarized 
in Table |3[ On the basis of xPT and preliminary unquenched resultsc2l, Sharpe estimates that S'C/(3)- 
breaking corrections | and unquenching may lead to an O (10%) increase in Bx and ascribes an 
O (15%) error to Bx to account for the uncertainties in this estimatelll. Bijnens et a/. t3 reach similar 
conclusions. I choose to include these effects as a contribution to the error but not to the central 
value. 

On the basis of these conclusions and the results given in Table |^, I quote: 
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nlo 



(2 GeV) = 0.65(10) or S^i° = 0.89(14) 



(13) 



where 5"^'° 



^ is the two-loop RG-invariant S-parameter obtained from Bj^ (2 G.eV) with nj = 3 and 
as(2GeV) = 0.3. Again, this result compares very favorably to the SM fit ofli3, but this time with 

as a fit parameter instead of fs^ ■' 
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r, 07+0.34 



"Calculations are performed with degenerate or nearly degenerate s and d quarks. 
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Table 4: Parameters of our calculations. "^7^ cfs" is the number of gauge-field configurations on which the various matrix 

elements are computed (i.e. our statistics), csw is the mean-field-improved coefficient of the SW term. a~^{mp) is the 

inverse lattice spacing as determined from a calculation of the p-meson mass. 

^P lattice size # cfs csw a~^(rrap) 

"6^2 24^ X 48 188 1.442 2.57(8) GeV 

6.0 16^ X 48 498 1.479 1.96(5) GeV 



4 Conclusions 

The lattice provides a means for calculating AB = 2 and AS = 2 matrix elements from first prin- 
ciples. A reliable determination of these matrix elements will be crucial for testing the SM with the 
forthcoming experiments on CP violating B decays. Moreover, the hadronic matrix elements which 
appear in supersymmetric extensions of the SM can also be considered Ej. In the next few years, more 
and more unquenched calculations will be performed, enabling a better quantification of quenching 
effects and eventually yielding fully unquenched results. 
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Appendix: details of the calculations 

We describe quarks with a mean- field- improved, Sheikholeslami-Wohlert (MFl SW) action. Compared 
to the standard Wilson action, the leading discretization errors are formally reduced by a factor of 
order as{a), and the mean- field-improvement may give additional numerical suppression. We perform 
calculations at two values of the cutoff, corresponding to couplings /? = 3/(27ras(a)) = 6.0 (coarser 
lattice) and 6.2 (finer lattice). This enables us to quantify discretization errors. The parameters of the 
simulations are summarized in Table ^. Note that our simulations have high statistics. Unfortunately, 
because of the very high numerical cost of including the feedback of quarks on the gauge-fields, both 
calculations are performed in the quenched approximation. 

Because a physical pion would feel the boundaries of the box in which we work and because the 
algorithms we use slow down rapidly for lighter quarks, we are restricted to work with quarks with 
masses on the order of ms/2 or more. Thus, to obtain results at the physical values of the u, d and 
s masses, we perform all calculations for three values of the light-quark mass roughly in the range 
ms/2 — > rris and extrapolate or interpolate to the physical mass values. 

Furthermore, the graininess of our lattice forbids us from working with quarks whose masses are 
much larger than m-c- Thus, we perform all calculations for five values of the heavy-quark mass around 
that of the charm and extrapolate to the physical 6-quark mass. 

Finally, because Wilson fermions break chiral symmetry explicitely, the left-left operator, Oa_f=2; 
mixes with four-quark operators of different chirality. We subtract these wrong chirality contributions 
and match to the MS-NDR scheme at one loop. To estimate the systematic error associated with our 
procedure, we vary the matching scale in the range 1/a — > it /a. For the S-parameters discussed here, 
dependence on this scale is very small. 
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